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ABSTRACT 

Context. Galaxy evolution scenarios predict that the feedback of star formation and nuclear activity (AGN) can drive the transforma¬ 
tion of gas-rich spiral mergers into (ultra) luminous infrared galaxies and, eventually, lead to the build-up of QS0/elliptical hosts. 
Aims. We study the role that star formation and AGN feedback have in launching and maintaining the molecular outflows in two 
starburst-dominated advanced mergers, NGC 1614 (Dl = 66 Mpc) and IRAS 17208-0014 (Dl = 181 Mpc), by analyzing the 
distribution and kinematics of their molecular gas reservoirs. Both galaxies present evidence of outflows in other phases of their ISM. 
Methods. We used the Plateau de Bure interferometer (PdBI) to image the CO(l-O) and CO(2-l) line emissions in NGC 1614 and 
IRAS 17208-0014, respectively, with high spatial resolution (0'.'5-l'.'2). The velocity flelds of the gas were analyzed and modeled to 
And the evidence of molecular outflows in these sources and characterize the mass, momentum, and energy of these components. 
Results. While most (> 95%) of the CO emission stems from spatially resolved (~ 2 - 3 kpc-diameter) rotating disks, we also detect 
in both mergers the emission from high-velocity line wings that extend up to +500-700 km s“^ well beyond the estimated virial 
range associated with rotation and turbulence. The kinematic major axis of the line-wing emission is tilted by ~ 90° in NGC 1614 
and by ~ 180° in IRAS 17208-0014 relative to the major axes of their respective rotating disks. These results can be explained by 
the existence of non-coplanar molecular outflows in both systems: the outflow axis is nearly perpendicular to the rotating disk in 
NGC 1614, but it is tilted relative to the angular momentum axis of the rotating disk in IRAS 17208-0014. 

Conclusions. In stark contrast to NGC 1614, where star formation alone can drive its molecular outflow, the mass, energy, and 
momentum budget requirements of the molecular outflow in IRAS 17208-0014 can be best accounted for by the existence of a so 
far undetected (hidden) AGN of Lagn ~ 7 x 10^^ L©. The geometry of the molecular outflow in IRAS 17208-0014 suggests that the 
outflow is launched by a non-coplanar disk that may be associated with a buried AGN in the western nucleus. 

Key words. Galaxies: individual: IRAS 17208-0014, NGC 1614 - Galaxies: ISM - Galaxies: kinematics and dynamics - Galaxies: 
nuclei - Galaxies: Seyfert - Radio lines: galaxies 


1. Introduction 

Different theoretical models and numerical simulations have 
long foreseen that massive gas outflows powered by active galac¬ 
tic nuclei (AGN) and/or nuclear starbursts could link the growth 
of black holes and their hosts during the history of galaxy 
evolution, as inspired by the ground-breaking work of Silk & 
Rees (1998). Outflows could prevent galaxies from becoming 
overly massive and help regulate the fueling of both the star 


Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany), and IGN (Spain). 


formation and the nuclear activity (e.g.. King 2003; Di Mat- 
teo et al. 2005, 2008; Ciotti & Ostriker 2007; Narayanan et 
al. 2008; Silk & Nusser 2010; Ishibashi & Fabian 2012; Faucher- 
Giguere & Quataert 2012; Hopkins et al. 2012, 2014; Zubovas 
& King 2012, 2014; Gabor & Boumaud 2014). 

There is mounting observational evidence of gas outflows 
in different populations of starbursts and active galaxies, in¬ 
cluding ultra luminous infrared galaxies (ULIRGs), radio galax¬ 
ies, quasars, and Seyferts. The outflow phenomenon concerns 
virtually all the phases of the interstellar medium (ISM), in¬ 
cluding the ionized gas (e.g., Heckman et al. 1990; Colina et 
al. 1991; Martin 1998; Holt et al. 2008; Crenshaw et al. 2010; 
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Rupke & Veilleux 2013; Arribas et al. 2014), the neutral atomic 
medium (e.g., Heckman et al. 2000; Rupke et al. 2002, 2005a, 
2005b, 2005c; Morganti et al. 2005; Martin 2005, 2006; Chen 
et al. 2010), as well as the molecular gas phase (e.g., Fer- 
uglio et al. 2010; Fischer et al. 2010; Sturm et al. 2011; Alat- 
alo et al. 2011, 2015; Chung et al. 2011; Aalto et al. 2012; 
Dasyra & Combes 2012; Maiolino et al. 2012; Combes et 
al. 2013; Morganti et al. 2013; Veilleux et al. 2013; Cicone 
et al. 2012, 2014, 2015; Garcia-Burillo et al. 2014; Gonzalez- 
Alfonso et al. 2014; Davies et al. 2014; Emonts et al 2014; 
Sakamoto et al. 2014). This multi wavelength evidence under¬ 
lines that the overall impact of outflows can only be quantifled 
if all of their ISM phases are studied in different populations of 
galaxies. 

Although ionized gas flows are seen to be typically a few 
100 km s“^ faster than their atomic or molecular counterparts, 
the neutral medium is seen to evacuate mass (Mout) and energy 
(E’en!) at comparably higher rates {MQ^x{neutral) ~ (10 - 100) x 
Moyii(ionized) and EQ^djieutral) ~ (1 - 100) x EQ^i{ionized)\ e.g., 
Rupke & Veilleux 2013). Since the molecular component is the 
most massive phase of the neutral ISM in the central kiloparsec 
regions, the study of molecular outflows is therefore crucial to 
gauging their global impact on galaxy evolution. 

The high spatial resolution and sensitivity capabilities, al¬ 
lied with the enhanced bandwidth of current mm-interferometers 
have revealed the existence of cold molecular outflows at 
kiloparsec scales in a growing number of (U)LIRGs. Out¬ 
flows are detected by the low-level broad line wings (up to 
~ + 1000 km s“^) identified mainly, yet not exclusively, in CO 
lines (Feruglio et al. 2010, 2013; Cicone et al. 2012, 2014). Trac¬ 
ers more specific to the dense molecular gas, such as low and 
mid-J rotational lines of HCN and HCO^, have been used to 
study the outflow in the prototypical AGN-dominated ULIRG 
Mrk 231 (Aalto et al. 2012, 2015a). 

Cicone et al. (2014) present new CO observations of seven 
ULIRGs and report detecting molecular outflows in four out 
of seven of their targets. The outflows, extending on kiloparsec 
scales, present mass load rates (dM/dO of several 100 M© yr“^ 
Based on their compiled sample of 19 sources, they found that 
star formation can drive outflows with dM/dr values up to two 
to four times the value of the star formation rate {SER). How¬ 
ever, the presence of an AGN, even if weak (LACN/^boi <0.10), 
can strongly boost an outflow with dM/dr oc Lagn- In the most 
extreme quasars (LACN/^boi ^ 0.80), outflow rates can be up to 
100 X 5 ER, and therefore have a decisive impact on the fueling 
of the star formation and nuclear activity of their hosts. 

In an attempt to study the specific role that star formation 
versus AGN feedback have in launching and maintaining molec¬ 
ular outflows in mergers, we used the Plateau de Bure inter¬ 
ferometer (PdBI) to image the CO emission in two starburst- 
dominated (U)LIRGs, NGC 1614 and IRAS 17208-0014, with 
high spatial resolution (0'.'5-l'.'2). The two targets present pre¬ 
vious evidence of an outflow in several gas tracers, including 
the ionized, the atomic, and the molecular ISM (NGC 1614: 
Bellocchi et al. 2012; IRAS 17208-0014: Sturm et al. 2011; 
Rupke & Veilleux 2013; Arribas et al. 2014). Both systems are 
classified as mergers in an advanced stage of interaction, and 
they host massive episodes of star formation. Interestingly, al¬ 
though recently questioned by Aalto et al. (2015b) in the case 
of IRAS 17208-0014, there has been no evidence to date for an 
energetically relevant AGN in these mergers at their late stage 
of evolution (LACN/^boi < 0-1 “ 11%; e.g., Gonzalez-Martm 
et al. 2009; Rupke & Veilleux et al. 2013; Herrero-Illana et 
al. 2014). The absence of any clear signs of ongoing nuclear 


activity makes the choice of these galaxies optimum in order 
to test the canonical evolutionary scenario where gas-rich spiral 
mergers lead first to a (U)LIRG and, eventually, under the action 
of AGN feedback to the build-up of QS0/elliptical hosts (e.g., 
Sanders et al. 1988; Veilleux et al. 2009). Furthermore, the high 
resolution of the PdBI observations used in this work are instru¬ 
mental in constraining the geometry, as well as the mass, mo¬ 
mentum, and energy budgets of the molecular outflows unveiled 
in NGC 1614 and IRAS 17208-0014. 

We describe in Sect. 2 the two targets that are the subject of 
this study. Section 3 describes the PdBI observations. Section 4 
presents the CO line maps obtained and discusses the observa¬ 
tional evidence of molecular outflows in both sources. Section 5 
describes the basic properties of the molecular outflows. We dis¬ 
cuss these results in the context of the molecular outflow phe¬ 
nomenon in (U)LIRGs in Sect. 6. The main conclusions of this 
work are summarized in Sect. 7. 

2, Targets 

2.1. NGC 1614 

NGC 1614 (Da = 64 Mpe, Dl = 66 Mpe; 1" = 310 pc) is an 
Hll-classifled LIRG (Ljr = 4.5x 10^^ L©) in an advanced stage of 
interaction. The bolometric luminosity of NGC 1614 is mainly 
powered by an extreme star formation episode with an integrated 
S ER of ~ 50M©yr“^ (Alonso-Herrero et al. 2001, U et al. 2012). 
Star formation is fed by a ~ 3 x 10^ M© molecular gas reser¬ 
voir (Konig et al. 2013). Numerical simulations suggest that the 
burst was triggered after a minor merger (Vaisainen et al. 2012). 
In their multi wavelength study, Herrero-Illana et al. (2014) con¬ 
cluded that there is no evidence of an AGN in the nucleus of 
NGC 1614 down to an AGN-to-bolometric luminosity ratio of 
< 10%. In agreement with this picture, not detecting the nuclear 
source in the 435yum continuum image of the galaxy obtained 
with ALMA led Xu et al. (2015) to rule out a Compton-thick 
AGN in NGC 1614. 

The molecular gas in this source, imaged by the OVRO, 
SMA, CARMA, and ALMA interferometers, is distributed in a 
nuclear disk and a fainter extension (Wilson et al. 2008; Olsson 
et al. 2010; Konig et al. 2013; Sliwa et al. 2014; Xu et al. 2015; 
Usero et al. in prep.). The recent 0'.'5-CO(2-1) map of the SMA, 
published by Konig et al. (2013), spatially resolved the nuclear 
disk into an asymmetric ring of ~ 230 pc radius, similar to the 
CO(6-5) image of the central region obtained with ALMA by 
Xu et al. (2015). The disk shows an extension to the northwest 
in the SMA image. However, the lower resolution (3") CO(l-O) 
map of CARMA (Sliwa et al. 2014) shows an extension to the 
east. 

The HST-NICMOS observations published by Alonso- 
Herrero et al. (2001) show a compact starburst nucleus (with an 
estimated age of > 10 Myr) surrounded by a ~ 2" (600 pc) diam¬ 
eter ring of supergiant HII regions revealed by Patr line emission, 
which is associated with a younger starburst (5-10 Myr). Lower 
level emission is also detected farther out to the east over a spiral 
structure out to r ~ 10" (~ 3 kpc). The luminosities of the giant 
HII regions in the ring are an order of magnitude brighter than 
30Dor. The nuclear ring shows an extremely high value of the 
Ha surface brightness: ~ 60 x 10^^erg s“^kpc“^ or, equivalently, 
a SER surface density of ~ 30 M©yr“^kpc"^. Such values are 
comparable to those derived across several kpc in high-z sub¬ 
millimeter galaxies (Menendez-Delmestre et al. 2013). The star¬ 
forming ring is also prominent in the high-resolution radio con¬ 
tinuum images (Olsson et al. 2010; Herrero-Illana et al. 2014) 
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and in the mid-infrared PAH emission maps of the galaxy (Diaz- 
Santos et al. 2008; Vaisainen et al. 2012). 

The star formation properties of the molecular gas in 
NGC 1614 stand out as non-standard on several counts. Among 
the objects studied by Garcia-Burillo et al. (2012), NGC 1614 
shows the most extreme (high) star formation efficiency of the 
dense molecular gas traced by the 1-0 line of HCN (S T’^’dense = 
Tfir/T^cn “ ^o/(K km s“^ pc^)). Furthermore, Xu et 

al. (2015) find a breakdown of the Kennicutt-Schmidt (KS) law 
on the linear scale of 100 pc in the nuclear ring, based on the gas 
column densities derived from the CO(6-5) map. 

The kinematics of the ionized gas traced by the Her line, 
which was studied by Bellocchi et al. (2012), show significant 
departures from circular motions. Bellocchi et al. (2012) find a 
main component associated with a rotating gas disk and a second 
(broader) component blueshifted by up to 300 km s“^ relative to 
the systemic velocity (fsys) of the galaxy. The latter was inter¬ 
preted as the signature of a dusty outfiow in the ionized gas. The 
blueshifted emission then comes from the outfiow lobe oriented 
toward us, and the far-side lobe remains undetected at optical 
wavelengths owing to extinction. The projections of the kine¬ 
matic major axes show that the rotating disk and the outfiow 
axis differ by ~ 90°, an indication that the outfiow is oriented 
along the minor axis of the disk. The estimated outfiow rate for 
the ionized gas is ~ 44Mo yr“^ (Colina, private communication). 

2.2. IRAS 17208-0014 

IRAS 17208-0014 (Da = 167 Mpe, Dl = 181 Mpe; 1"= 
810 pc) is a coalesced merger classified as an HII ULIRG 
(Lir = 2.4 X 10^^ Lq) with the multiwavelength characteris¬ 
tics of an obscured nuclear starburst (SFR ~ 240 M© yr“^) as 
the main driving source of its bolometric luminosity (Momjian 
et al. 2003; Nardini et al. 2009, 2010; Rupke & Veilleux 2013; 
Teng & Veilleux 2010; Iwasawa et al. 2011). However, the mor¬ 
phology of the X-ray source, analyzed by Iwasawa et al. (2011), 
shows hard X-ray emission and the tentative detection of a 
high-ionization Fe Ka line toward the NED nuclear position. 
Gonzalez-Martm et al. (2009) classify the source as a candi¬ 
date Compton-thick AGN with Lx ~ 10"^^ erg s“\ leaving room 
for an obscured but not energetically dominant AGN in the nu¬ 
cleus of IRAS 17208-0014: LACN/^boi ~ 0.1%. However, Rupke 
& Veilleux (2013) estimate this fraction to be ~ 11% based on 
MIR diagnostics. The MIR 8.7yum continuum emission appears 
clearly extended over ~ 1.6 kpc. The galaxy also shows extended 
11 3Jim PAH emission, but the equivalent width of this PAH fea¬ 
ture indicates the presence of an additional nuclear continuum 
source that might be associated with the obscured AGN (Alonso- 
Herrero et al. 2014). 

The optical morphology of the galaxy is patchy because of 
dust obscuration, but there is a well-defined nuclear source em¬ 
bedded in an extended disk of ~2 kpc-radius detected at near in¬ 
frared (NIR) wavelengths (Scoville et al. 2000). The NIR disk 
has a elliptical-like profile. The outer disk of the galaxy 
shows two tidal tails, which betray the past interaction (Mel- 
nick & Mirabel 1990; Murphy et al. 1996). The more recent 
high-resolution AO-assisted NIR observations done by Keck 
(Medling et al. 2014) have resolved the inner disk into two over¬ 
lapping nuclear disks with a small separation of ~ 072 (200 pc). 

Star formation is currently fueled by a massive molecular gas 
reservoir (~ (6 - 7) x 10^ M©: Downes & Solomon 1998; Wil¬ 
son et al. 2008). The molecular disk of IRAS 17208-0014 was 
imaged in the CO(l-O) line with the PdBI with a ~ 3'.'5-spatial 
resolution by Downes & Solomon (1998). The 2" (1.6 kpe)- 


diameter disk, barely resolved in these observations, shows a 
rotating pattern. The ~ l"-CO(3-2) SMA map of Wilson et 
al. (2008) reveals a compact ~ 175 (1.2 kpc) diameter disk with 
fainter extensions farther out to the northwest and to the east. 

Superposed on the overall rotating pattern, which dominates 
the kinematics of the disk, there is multiwavelength evidence 
of an outfiow found in the ionized and in the neutral gas. In 
spite of previous negative evidence of an outfiow in the ion¬ 
ized gas (Westmoquette et al. 2012; Rupke & Veilleux 2013), 
Arribas et al. (2014) report the detection of an outfiowing com¬ 
ponent, based on very high S/N Ha spectroscopy. The maximum 
velocities of the ionized outfiow reach moderate values of up 
to i; - Fsys ~ -300 km s“^ and the estimated outfiow rate for 
the ionized gas is ~ 16 M© yr"^ (Colina, private communica¬ 
tion). There is additional evidence of an outfiow in the neutral 
gas in IRAS 17208-0014. Flows with maximum velocities of 
up to F - t>sys ~ -650 kms“^ and mass rates dM/dt ~ 34 M© 
yr"^ were identified in the atomic gas in a spatially extended 
(r ~ 1.6 kpc) structure traced by Nal D absorption lines (Rupke 
& Veilleux 2013). Furthermore, outfiow velocities up to v-v^ys ~ 
-370 km s“^ with associated outfiow rates of ~ 90 M© yr“^ were 
also identified in the cold molecular gas through OH absorption 
lines (Sturm et al. 2011). However, this absorption line study, 
done with a spatial resolution of ~ 10", was not optimized to 
derive the size and the geometry of the molecular outfiow struc¬ 
ture, leaving the true outfiow rate largely unconstrained. More 
recently, the optical and near-infrared emission line ratios ana¬ 
lyzed by Medling et al. (2015) have confirmed that there is a 
multiphase outfiow of shocked gas extending over a wide range 
of spatial scales from r ~ 400 pc to r ~ 5 kpc in IRAS 17208- 
0014. 

3. Observations 

3.1. NGC 1614 

Observations of NGC 1614 were carried out with the PdBI array 
(Guilloteau et al. 1992) between June 2012 and March 2013. We 
used the ABD configurations and six antennae. We observed the 
/=l-0 line of CO (115.271 GHz at rest). During the observations 
we used the two polarizations of the receiver and the 3.6GHz- 
wide WideX spectral correlator of the PdBI; this is equiva¬ 
lent to 9500 km s“^ at the working frequency. Rest frequencies 
were corrected for the recession velocity initially assumed to be 
Vo(HEL) = 4795 km s“^ The systemic velocity, redetermined 
in this work, is v^y^(HEL) = 4763 km s“^ Observations were 
conducted in single-pointing mode with a field of view (primary 
beam size) of 43" centered at cr 2 ooo = 04^34^^00.03^ and ^2000 = 
-08°34'44.57". The latter is very close to the nominal positions 
of the optical and near infrared centers, as determined by Neff 
et al. (1990) and Skrutskie et al. (2006). We redetermined the 
dynamical center of the galaxy in this work, which coincides 
within the errors with the central peak detected in the 5 GHz 
radio continuum images obtained by Olsson et al. (2010) with 
MERLIN ( 0^2000 = 04^34'”00.03^ and ^2000 = -08°34'45.01"). 
Visibilities were obtained through on-source integration times of 
20 minutes framed by short (~2 min) phase and amplitude cali¬ 
brations on a nearby quasar. The absolute flux scale in our maps 
was derived to 10% accuracy based on the observations of pri¬ 
mary calibrators whose fluxes were determined from a combined 
set of measurements obtained at the 30m telescope and the PdBI 
array. The bandpass calibration is accurate to better than 5%. 

The image reconstruction was done with the standard 
IRAM/GILDAS software (Guilloteau & Lucas 2000). We used 
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both natural (NA) and uniform (UN) weighting to generate two 
different versions of the CO line map with a size of 51" and 
OV 1/pixel sampling. The corresponding synthesized beams are 
3'.'8 X 1'.'2 at PA =12° for the NA weighting map, and 1'.'5 x O'.'S 
at PA = 23° for the UN weighting map. The conversion fac¬ 
tors between flux (Jybeam“^) and temperature units (K) are 
20 K Jy“^ beam and 82 K Jy“^ beam for the NA and UN weight¬ 
ing datasets, respectively. The corresponding point source sensi¬ 
tivities derived from emission-free channels are 0.8 mJy beam“^ 
(NA) and 0.9 mJy beam“^ (UN) in 9.8 MHz (26 km s“^)-wide 
channels. The actual noise levels increase by up to a factor 2-3 
in channels with strong line-emission due to inherent limitations 
in the deconvolution process. An image of the continuum emis¬ 
sion of the galaxy was obtained by averaging those channels free 
of line emission. The point source sensitivity for the continuum 
map is 0.1 mJy beam“^ in the highest resolution version of the 
data obtained by adopting UN weighting. 

3.2. IRAS 17208-0014 

Observations of IRAS 17208-0014 were carried out with the 
PdBI in March 2012. We used the B configuration and six 
antennae of the array. We observed the /=2-l line of CO 
(230.538 GHz at rest). As for NGC 1614, we used the two po¬ 
larizations of the receiver and the WideX correlator, and this 
is equivalent to 4900 km s"^ at the CO(2-l) frequency. Rest 
frequencies were corrected for an adopted recession velocity 
of Vo(HEL) = 12838 km s“^ We derived the systemic veloc¬ 
ity in this work, which is ~30 km s“^ blueshifted relative to 
Vo'. v^y^iHEL) = 12808km s“^ We used a single-pointed ob¬ 
serving mode with a primary beam size of 23" centered at 
0^2000 = 17^23'^21.90^ and ^2000 = -00°17'00.90", which is 
< 1" away from the radiocontiuum peak detected in the 1.4 GHz 
VLB I maps of Momjian et al. (2003) and the dynamical cen¬ 
ter of the galaxy, which was determined in this work ( 0^2000 = 
17^23'”21.96" and ^2000 = -00°17'00.87"). Visibilities were ob¬ 
tained through on-source integration times of 20 minutes framed 
by short (~2 min) phase and amplitude calibrations on a nearby 
quasar. The absolute flux scale in our maps was derived to 10% 
accuracy, and the bandpass calibration is accurate to better than 
5%. 

We performed image deconvolution using IRAM/GILDAS 
software. As for NGC 1614, we used both NA and UN weight¬ 
ing to generate two versions of the CO map with a size of 26" 
and OV 1/pixel sampling. The corresponding synthesized beams 
are 078 x 075 at RA = 18° for the NA weighting map, and 
076x075 at RA = 61° for the UN weighting map. The conversion 
factors between flux (Jybeam“^) and temperature units (K) are 
62 K Jy“^ beam and 84 K Jy"^ beam for the NA and UN weight¬ 
ing datasets, respectively. The corresponding point source sensi¬ 
tivities derived from emission-free channels are 0.9 mJy beam“^ 
(NA) and 1.2 mJy beam“^ (UN) in 20 MHz (27 km s“^)-wide 
channels. The noise figures are to be increased by a factor 2-3 
in channels with strong line emission. A continuum image of the 
galaxy was obtained by averaging channels free of line emission 
with a point-source sensitivity of 0.1 mJy beam“^ in the map de¬ 
rived with UN weighting. 

4. Results 

4.1. NGC 1614: channel maps and emission centroids 

Figure 1 shows the CO(l-O) velocity-channel maps observed in 
NGC 1614 with a spatial resolution of 378 x 172 at PA = 12°, 


obtained using natural weighting of the data to recover the max¬ 
imum percentage of low-level emission in extended structures. 
We display channel maps grouped by pairs with a velocity spac¬ 
ing of 60 km s“^ from v - fsys = 0 km s“^ with i;sys(HEL) = 
4763 km s“^ a value determined by the best-fit kinemetry solu¬ 
tion found for the mean velocity field described in Sect. 4.2. 

Virtually all the CO line emission detected over significant 
3(T-levels comes from the central r ~ 10" (3 kpc) region that 
is displayed in Fig. 1. Overall, this emission extends over a 
roughly symmetric ~ 900 km s“^-wide velocity span: v - t>sys ~ 
[-1-450, -450] km s“^ However, the bulk of the emission, up to 
^ 98% of the grand total, comes from a more restricted velocity 
range: | v - t>sys l< 210 km s“\ hereafter referred to as the line 
core. 

Figure 2 shows the CO spectrum spatially integrated over 
a square region of 4 kpc size to cover the full extent of the 
NGC 1614 disk. The emission from the line core can be sat¬ 
isfactorily fit by three Gaussian components. In this fit we as¬ 
signed three initial values for the velocity centroids of the Gaus¬ 
sian components: one around t>sys. and two velocities symmetri¬ 
cally offset by +100 km s“^ with respect to t>sys to cover the ex¬ 
pected symmetric range caused by rotation found in the model 
of Bellocchi et al. (2012). However, outside the line core there 
is an excess of fainter but nevertheless statistically significant 
emission (up to ~ Her in velocity-integrated units) which ex¬ 
tends up to I - t>sys l~ 450 km s“\ hereafter referred to as the 
line wings (i.e., 210 km s“^<| v - t>sys 1^ 450 km s“^). Emission 
stemming from the line core shows the signature of a rotating 
disk with a kinematic major axis oriented roughly north-south: 
emission at redshifted (blueshifted) velocities comes from the 
northern (southern) side of the disk, as can be guessed by in¬ 
specting Fig. 1. By contrast, outside the line core, the kinematic 
major axis for the line wings seems to tilt progressively toward 
the east-west axis (see Fig. 1). 

In an attempt to quantify the kinematic decoupling of the 
emission of the line wings relative to the line core, we have de¬ 
rived the centroids of CO emission in NGC 1614 in [Aa, AS]- 
units. The left-hand panel of Fig. 3 shows the CO centroids for 
the same velocity channels displayed in Fig. 1. The centroids 
of CO(l-O) emission are derived using as weighting function a 
power law of the flux measured at a given offset within the inner 
20" field of view with an index n = 1 /2 in order to optimize 
the contribution from extended emission and reduce accordingly 
the noise caused by dumpiness in this determination. Emission 
centroids calculated for velocity channels which are equidistant 
from fsys are connected so as to visualize the change in PA as 
a function of velocity, shown in the middle panel of Fig. 3. The 
PA is defined for the redshifted component of the connected cen¬ 
troids, and it is measured east from north (i.e., negative values to 
the west). The right-hand panel of Fig. 3 illustrates the change 
in average radius r, defined as half of the mutual distance be¬ 
tween the connected centroids, as a function of velocity. The 
left-hand panel of Fig. 3 identifies the dynamical center at [Aa, 
A^]=[0", 0"] as determined in Sect. 4.2. The centroid barycenter 
is slightly shifted relative to the dynamical center because of the 
non-axisymmetries in the brightness distribution of CO. 

As illustrated in the middle panel of Fig. 3, the value derived 
for PA changes by ~ 80° from | t>-t>sys I = 60 km s“\ where PA ~ 

+ 13°,to| v-v^ys I = 360 km s“^ where FA - 64°. We estimate 

a weighted mean for PA in the line core of FAcore = 2 + 5°, and 
this is close to the value derived for FA using the kinemetry fit 

of Sect. 4.3 (FA-8° + 5°). The corresponding weighted mean 

for FA in the line wing region is FA wing = -68 + 15°. We can 
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Fig. 1. CO(l-O) velocity-channel maps observed with the PdBI in the nucleus of NGC 1614 with a spatial resolution of 378 x 172 at PA = 12° 
obtained with natural weighting (beam is plotted as a fllled ellipse in the bottom left corner of the lower left panel). Velocity resolution is 
60 km s“^ We show a held of view of 20", i.e., ~l/2 the diameter of the primary beam at 113.5 GHz. We display channel maps grouped by pairs 
from V - Vsys - 0 km s~^ in steps of 60 km s~^ with i>sys(HEL) = 4763 km s“^ from -480 km s“^ to 480 km s“^. Emission at blueshifted velocities is 
displayed in gray scale and black contours, and emission at redshifted velocities is displayed in red contours. Contour levels are 3cr, 5cr, 8cr, 12cr, 
and 20cr with Icr = 0.5 mJy beam"^ for the channels of the line wings (|i> - i>sysl > 210 km s"^), and 3cr, 5cr, 8cr, 12cr, 20cr to 80cr in steps of 20cr, 
with Icr = 1.5 mJy beam"^ for the channels of the line core (}v - i;sysl < 210 km s“^; boxes highlighted in blue). The position of the dynamical 
center ([Acr, A 6 ] = [0",0"] = [cr 2 ooo = 04^34'”00.03^ 62000 = -08°34'45.01"]) is highlighted by the (blue) square marker. 


therefore conclude that the evidence of a kinematic decoupling, 

reflected in the measured difference Z^A^ing - P^cotq -70 +16° 

is a significant ^ 4Acr result. However, the right-hand panel of 
Fig. 3 shows that there is no statistically significant change in 
the radial distance, which stays fairly constant at a value of r ~ 
175 - 2" (0.5 - 0.6 kpc) as a function of the velocity offset. 


4.2. NGC 1614: moment maps 

Figure 4 shows the first three moment maps derived from the 
CO(l-O) line data obtained by applying uniform weighting to 
the visibilities inside the velocity interval v - L>sys = [+600,- 
600] km s“^. Uniform weighting allows us to reach a spatial reso¬ 


lution of 175x078 at PA = 23°. The velocity-integrated intensity 
map shown in the left-hand panel of Fig. 4 was derived assuming 
a 3(T-clipping on the intensities. The mean-velocity field and the 
velocity width maps, shown in the middle and right-hand panels 
of Fig. 4, respectively, were both derived by assuming a more 
demanding 5(T-clipping on the intensities to maximize the relia¬ 
bility of the images. 

The left-hand panel of Fig. 4 shows that the CO emission 
comes from a spatially-resolved 11" x 8" (3.4 kpc x 2.5 kpc) 
lopsided molecular disk. The disk shows an arc-like feature that 
runs from the northeast to the southeastern side. The strongest 
emission peak is located ~ 300 pc northwest relative to the cen¬ 
tral (0", 0")-offset. 
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Fig. 2. CO(l-O) spectrum obtained by spatially integrating the emission over the disk of NGC 1614. We show the three individual Gaussian 
components fitted to the spectrum (gray curves) and the implied overall fit combining them (red curve) (left panel). A close-up view of the overall 
fit to the spectrum is shown in the middle panel to highlight emission from the line wings. Residuals and associated Gaussian fits are shown in the 
right panel with the typical Icr level indicated by the dotted line. 


NGC1614 




Fig. 3. Centroids of CO(l-O) emission in NGC 1614 (in [Acr, A(5]-units) derived with a velocity resolution of 60 km s ^ using a power law of 
the flux as weighting function with index n - 1/2, are displayed in the left panel . Emission centroids calculated for velocity channels that are 
equidistant from i>sys (i.e., with a common \v - i>sys|-offset) are connected to visualize the change of PA, defined for the redshifted component of 
the connected centroids and measured east from north, as a function of velocity {middle panel). We also plot the change in the average radius r, 
defined as half of the mutual distance from the connected centroids, as a function of velocity {right panel). Red color is used for velocity channels 
of the line wing and black color for the line core. The dotted black (red) line in the middle panel identifies the + Icr range around the average PA 
of the disk as derived inside the line core (line wing): PAcore = 2 + 5° (PA^mg = -68 + 15°). The green marker in the left panel identifies the 
dynamical center. Error bars as a function of velocity channel are shown in the middle and right panels. Error bars in the left panel are for the line 
core. 


The middle panel of Fig. 4 confirms the signature of a spa¬ 
tially resolved rotating disk, as described in Sect. 4.1. When 
assuming that the spiral features seen in the optical and NIR 
images of the galaxy are trailing, the gas rotation should be 
counterclockwise to account for the observed velocity gradi¬ 
ent. This constrains the geometry of the galaxy disk with the 
eastern (western) side being the far (near) side. The choice of 
the (0", 0")-olFset as the nucleus maximizes the symmetry of 
the mean-velocity field. This location corresponds to the maxi¬ 
mum velocity witdh (FWHM) measured in the disk, shown in the 
right-hand panel of Fig. 4, and also to the radio continuum peak 
of 01sson et al. (2010). Therefore this position qualifies as the 
dynamical center of the galaxy. The molecular gas is very turbu¬ 
lent: we estimate an average FWHM-80 km s“^ over the disk. 
The value of FWHM is lowered to -70 km s“^ if we exclude the 


central r = 1.5" (0.5 kpc) of the disk to minimize the effects of 
beam smearing for our current spatial resolution. This implies an 
intrinsic velocity dispersion of cr ^ 30 km s“^ on average. 

An inspection of the middle panel of Fig. 4 indicates that 
while the kinematic major axis is oriented roughly north-south, 
the kinematic minor axis is not oriented east-west, but rather at 
an angle (> 15°). This global tilt proves that there are signifi¬ 
cant departures from circular rotation. Additional wiggles in the 
velocity field, identified on the eastern side of the disk, betray 
non-circular motions on smaller spatial scales. 

The left-hand and middle panels of Fig. 5 overlay the CO line 
intensities on the HST NICMOS images of NGC 1614 obtained 
in the Patr line and in the 1 .6pm continuum by Alonso-Herrero et 
al. (2001). While a sizable fraction of the CO emission originates 
in the nuclear star-forming ring, which is barely resolved by the 
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Fig. 4. Left panel. CO(l-O) integrated intensity map obtained in the disk of NGC 1614. The map is shown on color scale with contour levels 
5%, 10%, 15%, 20%, to 90% in steps of 10% of the peak value =14 Jy km s"^beam"^ (~ 50cr). Middle panel The CO(l-O) isovelocity contours 
spanning the range (-120 km s“^ 120 km s"^) in steps of 20 km s"^ are overlaid on a false-color velocity map (linear color scale as shown). 
Velocities refer to i;sys(HEL) = 4763 km s“^ Right panel Overlay of the CO(l-O) line widths (FWHM) shown in contours (30, 45 to 195 km s"^ 
in steps of 30 km s"4 on a false-color width map (linear scale as shown). The fllled ellipses in the bottom left comer of each panel represent the 
CO(l-O) beam size obtained with uniform weighting (1'.'5 x O'.'S at PA - 23°). The position of the dynamical center is highlighted by the cross 
marker. 
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Fig. 5. Left panel Overlay of the CO(l-O) intensity contours (levels as in the left-hand panel of Fig. 4) on the Pacr emission HST map of NGC 1614 
(color scale as shown in arbitrary units). Middle panel Same as left panel but with CO contours overlaid on the HST/NICMOS l.6pm continuum 
image (color scale as shown in arbitrary units). Right panel Same as middle panel but zooming in on the continuum emission image of NGC 1614 
at 113.5 GHz, in contours (10% to 90% in steps of 10% of the peak value = 2.4 mJy beam"^ ~ 24cr), overlaid on the HST NICMOS image. The 
fllled ellipses represent the CO beam size. 


PdBI, fainter CO emission is also detected extending farther out 
to the east up to r ~ 6" (1.9 kpc), where strong Pacr emission 
identifies the star-forming complexes of the eastern spiral arm. 
The right-hand panel of Fig. 5 shows that, unlike the CO line 
emission, the mm-continuum emission at 113.5 GHz is restricted 
to the nuclear ring. 


4.3. NGC 1614: kinematic modeling 

4.3.1 . Kinemetry fit: the coplanar solution 

To quantify the deviations from circular motions discussed in 
Sect. 4.2, we used the mean-velocity field of Fig. 4 to find the 
best-fit Fourier decomposition of the line-of-sight velocities us¬ 
ing the software package kinemetry (Krajnovic et al. 2006). A 
key implicit assumption of this decomposition, which we ques¬ 


tion in Sect. 4.3.2, is that the gas kinematics in NGC 1614 can 
be modeled by orbits that lie at all radii in a coplanar geometry. 

Figure 6, which shows the radial profiles of the main terms 
of the Fourier decomposition of the velocity field, summarizes 
the results of this analysis. The ci term, which accounts for 
the (projected) circular component of the velocity field (t>circ )5 
is the dominant contributor to the observed motions, as illus¬ 
trated in Fig. 6a. Figure 7 shows the rotating disk solution of 
kinemetry and the associated residuals. The geometry of the ro¬ 
tating disk is defined by a position angle PA = -8 + 5 = 352 + 5° 
(measured east from north for the receding side of the ma¬ 
jor axis) and an inclination i = 33 + 2°. We also determined 
fsysCHEL) = 4763 + 10 km s“\ as part the best-fit solution. Rel¬ 
ative to Vcirc, non-circular motions (fnonc) derived up to the third 
order of the Fourier decomposition, are within a significant 10% 
to 30% range, as shown in Fig. 6b. Furthermore, Fig. 6c shows 
that the sign of the si term, which represents the (projected) ax- 
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Fig. 6. a) Cl radial profile (black curve) derived from the Fourier decomposition of the velocity field of NGC 1614. The ci term accounts for the 
(projected) circular component of the velocity field (fcirc)- We also plot (red curve) the radial profile of the (projected) non-circular motions (t^nonc) 
derived up to the third order of the decomposition, b) The radial profile of the t^nonc/^^circ ratio, c) The radial variation of the si/ci. The term 
represents the (projected) axisymmetric radial motions of the velocity field. Errorbars in all panels represent Icr errors. 
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Fig. 7. Left panel. Overlay of the velocity field of the best-fit rotating disk model (in contours spanning the range -100 km s"^ to 100 km s"^ in 
steps of 20 km s"^ on the CO(l-O) isovelocities (color scale as shown) of NGC 1614. Right panel Overlay of the CO(l-O) intensity contours 
(levels as in the left-hand panel of Fig. 4) on the residual mean-velocity field (in color scale as shown) obtained after subtraction of the best-fit 
rotating disk model from the observations. 


isymmetric radial motions of the velocity field, is negative up 
to r ~ 2" (0.6 kpc). This is indicative of inflowing motions at 
these radii for the adopted geometry of the disk when assuming 
the gas is coplanar. Deprojected, these in-plane motions would 
translate into a maximum amplitude of the radially inward fleld 
of <15 km s“^. These low amplitude inward motions can be ex¬ 
plained by density-wave-driven gas flows. 

4.3.2. A non-coplanar outflow solution 

While the coplanar solution discussed in Sect. 4.3.1 could ac¬ 
count for the observed trend of the mean velocities fitted by 
kinemetry, it nevertheless fails to explain the order of mag¬ 
nitude and, also, the 2D pattern of the high velocity emission 
found in the line wings, as argued in the following. 

The left-hand panel of Fig. 8 shows the position-velocity 
(p-v) plot taken along the kinematic major axis of the rotating 
disk {PA = 352°) with the best-fit (projected) rotation curve su¬ 
perposed. Most of the emission in the line core lies within the 


estimated virial range. The latter is defined by a combination 
of circular rotation (~ +100 km s“^), turbulence (FWHM/2~ 
±40 km s“^), and an estimated upper limit to the contribution 
from in-plane non-circular motions of ~ 1/2 x L>circ^ 
Colombo et al. 2014; Garcia-Burillo et al. 2014). However, emis¬ 
sion from the line wings, identified in Fig. 8, goes beyond the 
adopted virial range, arguing in favor either of a fast(er)-rotating 
disk or of a non-rotationally supported component. 

The middle panel of Fig. 8 shows the p-v plot taken along 
the kinematic minor axis {PA = 82°). The bulk of the emission 
from the line wings, which stems from the nuclear star-forming 
ring region, shows a distinct kinematic pattern: the high-velocity 
blueshifted (redshifted) gas is detected on the eastern (western) 
side of the ring up to r ~ 2" (0.6 kpc), as better shown in the 
close-up view of the right-hand panel of Fig. 8. Taken at face 
value, the observed pattern invalidates the interpretation of the 
line wing emission as due to a fast-rotating coplanar disk. A sim¬ 
ilar pattern can be traced continuously into the lower velocities 

^ This amounts to +50 km s"^ in NGC 1614. 
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Fig. 8. Left panel: Position-velocity (p-v) plot taken along the kinematic major axis {PA - 352°) of NGC 1614. The blue star markers delineate the 
(projected) rotation curve fitted in Sect. 4.3 and the allowed virial range around it (delimited by the blue lines). The red squares delimit the region 
where emission is detected from the line wings. Contour levels are -3cr, 3cr, 5cr, lOcr, 20cr, to 200cr in steps of 20cr with Icr = 0.8 mJy beam“k 
Middle panel: The p-v plot taken along the kinematic minor axis {PA - 82°). Contour levels are -3cr, 3cr, 5cr, lOcr, 20cr, to 140cr in steps of 20cr 
with Icr = 0.8 mJy beam"k The allowed virial range is delimited by the blue curves. Right panel: Same as middle panel but zooming in on the 
central +3" region around the center. The red squares delimit the region where emission is detected from the line wings. In all panels the velocity 
scales are relative to i>sys, and the spatial scales (X and Y) are relative to the dynamical center. We show the spatial resolutions projected along the 
major axis {left panel) and along the minor axis {middle and right panels) obtained with NA weighting. 


of the line core, which suggests that a single kinematic com¬ 
ponent is responsible for the velocity gradient observed along 
the minor axis. For the assumed geometry of the disk, this im¬ 
plies that the gas would be moving radially inward if it lay in the 
galaxy plane. In particular, emission from the line wings would 
be coming from gas moving inward at very high velocities of up 
to ~ (200 - 400)/vm(/) km s“^ = 370-740 km s~^. However, this 
regime is well beyond the expected range of inward velocities of 
a typical density-wave gas flow. 

In contrast, if we abandon the restriction of the coplanar so¬ 
lution for the gas, the kinematic pattern of the line wings can 
be satisfactorily explained either by a non-coplanar (less ex¬ 
treme) inflow solution or, alternatively, by a molecular outflow, 
depending on the adopted geometry. As argued below, the non- 
coplanar outflow solution is favored. In the outflow scenario, the 
blueshifted velocities would stem from the lobe pointing toward 
us on the eastern side of the disk, while the redshifted veloci¬ 
ties would be behind the disk on the western side. This outflow 
model for CO agrees perfectly with the geometry of the Ha out¬ 
flow discovered by Bellocchi et al. (2012). The non-detection 
of redshifted Ha emission indicates that this component is lo¬ 
cated behind the disk, and therefore remained undetected ow¬ 
ing to extinction. The extinction derived from the 1.6/im/2.2/im 
color map of Alonso-Herrero et al. (2001) appears to be higher 
on the western side of the ring. The similar spatial extent and 
(blueshifted) velocities measured for the Ha and CO features 
suggest that both are part of the same outflowing component. 
The observed tilt in the major axis of the molecular outflow 
relative to the PA of the rotating disk (~ 80° as discussed in 
Sect. 4.1) suggests that the outflow has a signiflcant component 
that is oriented perpendicular to the rotating disk. This geometry 
is typically found in star-forming-driven molecular outflows, like 
the ones discovered in the starbursts M 82 and NGC 253 (Garcia- 
Burillo et al. 2001; Walter et al. 2002; Bolatto et al. 2013b; Salas 
et al. 2014). However, we cannot exclude that an AGN can also 
contribute to drive the outflow (see discussion in Sect. 5.3). 


4.4. IRAS 17208-0014: channel maps and emission 

centroids 

Figure 9 shows the CO(2-l) velocity-channel maps observed 
in IRAS 17208-0014 with a spatial resolution of 076 x 075 at 
PA = 61°, obtained using uniform weighting of the data. We 
display channel maps grouped by pairs with a velocity spac¬ 
ing of 60 km s“^ from v - v^ys = 0 km s“^ with i;sys(HEL) = 
12808 km s“^ a value determined by the best-flt solution de¬ 
scribed in Sect. 4.6. 

Most of the CO emission stems from the central r ~ 2" 
(1.7 kpc) region shown in Fig. 9. The emission covers a wide 
and remarkably asymmetric velocity range: v - L>sys = [+450,- 
750] km s“^ However, up to ^ 99% of the total emission 
comes from a narrower and symmetric velocity range: v - 

~ [+450, -450] km s“^ Following the notation adopted in 
NGC 1614, we also refer to this central velocity interval as the 
line core and identify the emission outside this range as coming 
from the line wing (i.e., -450 km s“^> v - ^sys ^ -750 km s“^). 
Although at an admittedly lower level compared to the line core, 
emission from the line wing reaches an overall statistical signif¬ 
icance of about ~ 15(T in velocity-integrated units. 

Figure 10 shows the CO emission profile extracted from the 
position where the emission in the line wing is maximum at [Aa, 
A5] = [0",-0''2], which corresponds to ^2000 = 17'’23'"21.963", 
^2000 = -00°17'01.03". We used three Gaussian components 
to fit the emission from the line core. As in NGC 1614, the fit 
left a residual of emission, which in IRAS 17208-0014 is iden¬ 
tified at highly blueshifted velocities in the range v - L>sys = [- 
450,-750] km s“^ and a tentatively detected component around 
-950 km s“^ This blueshifted line wing has no redshifted coun¬ 
terpart, unlike in NGC 1614. Emission from the line core re¬ 
veals a rotating disk with a kinematic major axis oriented at a 
small angle relative to the east-west axis: emission at redshifted 
(blueshifted) velocities comes from the eastern (western) side of 
the disk, as illustrated in Fig. 9. Within the limits of our spa- 
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Fig. 9. C0(2-l) velocity-channel maps observed with the PdBI in the nucleus of IRAS 17208-0014 with a spatial resolution of 0'.'6 x 0'.'5 at 
PA = 61° (beam is plotted as a filled ellipse in the bottom left corner of the lower left panel). Velocity resolution is 60 km s"^ We show a field of 
view of 4", i.e., ~l/5 the diameter of the primary beam at 221.1 GHz. We display channel maps grouped by pairs from v - Vsys=0 km s"^ in steps of 
60 km s~^ with i;sys(HEL) = 12808 km s~^ from -720 km s"^ to 720 km s"^ Emission at blueshifted velocities is displayed in gray scale and black 
contours, and emission at redshifted velocities is displayed in red contours. Contour levels are 3cr, 5cr, and 8cr with Icr = 0.7 mJy beam"^ for the 
channels of the line wing (|i; - i>sysl > 450 km s"^), and 3cr, 5cr, 8cr, 12cr, 20cr, to lOOcr in steps of 20cr, and 160cr with Icr = 1.4 mJy beam"^ for the 
channels of the line core (|i> - Vsysl < 450 km s"^; boxes highlighted in blue). The position of the dynamical center ([Acr, A^] = [0",0"] = [cr 2 ooo = 
17^23"*21.96D 62000 = -00°17'00.87"]) is highlighted by the (blue) square marker. 
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Fig. 10. Same as Fig. 2 but for the CO(2-l) spectrum obtained toward the position that maximizes the emission in the line wing in IRAS 17208- 
0014 at 0^2000 = 17^23"*21.963D ^2000 = -00°17'01.03". 


tial resolution, Fig. 9 shows that emission in the line wing shifts 
progressively from the western to the eastern side of the disk 
as the velocity offset increases, i.e., showing an opposite trend 


compared to the line core. This is a signature of kinematic de¬ 
coupling of the gas emitting at high velocities. 

We have derived the centroids of CO emission in 
IRAS 17208-0014 to quantify the structure of the different kine- 
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Fig. 11. Same as Fig. 3 but for the centroids of CO(2-l) emission in IRAS 17208-0014. 


matic components following the same procedure as described in 
Sect. 4.1 for NGC 1614. Figure 11 shows both the centroids and 
the PA and r profiles obtained for the same velocity channels as 
are displayed in Fig. 9. As shown in the middle panel of Fig. 11, 
the value of PA stays fairly constant inside the line core around 
PAcore = 105 + 2°, an angle that is close to the value derived in 
Sect. 4.5 for the PA of the rotating disk (PA ~ 113° + 3°). In 
stark contrast, the PA values for the line wing centroids oscillate 
around a weighted mean of = 340° + 20°, an indication 

that its apparent kinematic major axis is virtually reversed: the 
PA of the line wing changes by an angle > 230° relative to the 
line core We conclude that the evidence of a kinematic de¬ 
coupling, refiected in the measured difference /^A^mg - PAcoyq ~ 
235 + 20° is a significant ^ lOcr result. Furthermore, the line 
wing channels show lower r values (~ O'.'l - 072 (80 - 160 pc)) 
compared to the line core (up to ~ 073 (240 pc)), indicating that 
the gas emitting at anomalous velocities lies at comparatively 
smaller radii. 

4.5. IRAS 17208-0014: moment maps 

Figure 12 shows the moment maps derived from the CO(2- 
1) line data inside the velocity interval v - v^ys = [-1000, 
500] km s“^. As for NGC 1614, we used different clippings: a 
3(T-clipping on the intensities to derive the zero-order moment 
(left-hand panel of Fig. 12) and a 5(T-clipping to derive the mean- 
velocity field and the velocity width maps (shown in the middle 
and right-hand panels of Fig. 12, respectively). 

The left-hand panel of Fig. 12 shows that the CO emission 
comes from a spatially-resolved 373 x 272 (2.7 kpc x 1.8 kpc) 
molecular disk feature that is noticeably elongated along PA ~ 
20°. The overall morphology of the CO disk is similar to the 
northeast-southwest elongated structure seen in the dust extinc¬ 
tion, as derived from the 2.2pm/l.Ipm HST color image of Scov- 
ille et al. (2000). However, compared to the CO(3-2) disk of Wil¬ 
son et al. (2008), the CO(2-l) disk is a factor 1.5-2 larger. As 
shown in the mean-velocity field of the middle panel of Fig. 12, 
which confirms the signature of a spatially resolved rotating 
disk, the morphological and kinematical major axes of the CO 
disk are roughly orthogonal. This implies that the intrinsic elon¬ 
gation of the disk should be more pronounced after deprojection. 
Three fainter protrusions and a detached clump extend the disk 

^ In the absence of a redshifted counterpart in the line wing, we derived 
the PA for this component using the central offset and the blueshifted 
channel centroids as reference. 


emission farther out up to r ~ 2.2" (1.8 kpc). Most of these 
CO extensions coincide in position with similar features of the 
2.2jim/l.lpm HST color image of Scoville et al. (2000). 

If, as discussed by Arribas & Colina (2003), the spiral arms 
and the most prominent tidal feature of the optical images of the 
galaxy (e.g., Melnick & Mirabel 1990; Solomon et al. 1997) are 
assumed to be trailing, gas rotation should be counterclockwise, 
and the observed CO velocity field would imply that the north¬ 
eastern side of the disk is the near side. The (0", 0")-offset, ini¬ 
tially adopted as the nucleus of the galaxy, maximizes the sym¬ 
metry of the mean-velocity field. It also corresponds to the maxi¬ 
mum velocity width (FWHM) measured in the disk, as shown in 
the right-hand panel of Fig. 12, and, also, to the peak of the con¬ 
tinuum emission at 221.1 GHz, as shown in the right-hand panel 
of Fig. 13. Therefore we adopt this position as the dynamical 
center of the galaxy. 

The average FWHM of the CO line in the disk is 
~170 km s“^ and this is lowered to ~150 km s“^ if we exclude 
the central r = 073 (0.2 kpc) of the disk to minimize the effects 
of beam smearing. This would imply an anomalously high av¬ 
erage velocity dispersion of cr ^ 64 km s“^ However, the same 
estimate is a factor of two lower if we derive it in the outer disk 
region (r > 175 (1.2 kpc)), an indication that there are unre¬ 
solved gradients of circular and non-circular motions in the disk. 
According to the middle panel of Fig. 12, the kinematic major 
axis of the rotating disk is oriented along PA ~ 110°. However, 
there are significant departures from circular rotation, identified 
by wiggles in the velocity field at several locations in the disk. 

The left-hand and middle panels of Fig. 13 overlay the 
CO line intensities on the HST optical and NIR images of 
IRAS 17208-0014 obtained in the F814W and F160W bands 
and published, respectively, by Arribas & Colina (2003) (from 
program 6346 by PI: K. Borne) and Scoville et al. (2000). The 
bulk of the CO emission is associated with the compact star- 
burst region, which hosts ^18 massive star clusters cataloged 
by Scoville et al. (2000). Unlike the CO line emission, the mm- 
continuum emission at 221.1 GHz, yet spatially resolved^, is 
mostly restricted to the compact starburst region, as shown in 
the right-hand panel of Fig. 13. 


^ The continuum emission has a deconvolved size of 073 x 072, deter¬ 
mined by task UVFIT of GILDAS. 


Article number, page 11 of 21 






















A&A proofs: manuscript no. sgarcia-burillo_26133_corrected 


IRAS17208-0014 


1 10 100 -200 -100 0 100 200 100 200 300 400 500 



□_^^^^^^_I_I_^_I_^^^_I_I_^^^^^^_I_I_d h_I_I_^^^^^^_I_I_^^^^^_I_I_^^^^^^_d d_^^_I_I_^^^^^^_I_I_^^^^^_I_I_^^^^_c 

2 1 0 , -1 -2 2 1 0 -1 -2 2 1 0 -1 -2 
Aa(arcsec) 

Fig. 12. Left panel. The CO(2-l) integrated intensity map obtained in the disk of IRAS 17208-0014. The map is shown in color scale with contour 
levels 1%, 2%, 5%, 10%, to 90% in steps of 10% of the peak value = 143 Jy km s"^beam"^ (-300 cr). Middle panel The CO(2-l) isovelocity 
contours spanning the range (-250 km s“^ 250 km s"^) in steps of 50 km s"^ are overlaid on a false-color velocity map (linear color scale as shown). 
Velocities refer to i>svs(HEL) = 12808 km s“^ Right panel Overlay of the CO(2-l) line widths (FWHM) shown in contours (50 to 550 km s"^ 
in steps of 50 km s"0 on a false-color width map (linear scale as shown). The filled ellipses in the bottom left comer of each panel represent the 
CO(2-l) beam size (0'.'6 x 0'.'5 at PA = 61°). The position of the dynamical center is highlighted by the cross marker. 
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Fig. 13. Left panel Overlay of the CO(2-l) intensity contours (levels as in the left-hand panel of Fig. 12) on the HST optical image of IRAS 17208- 
0014 obtained in the F814W band (color scale as shown in arbitrary units). Middle panel Same as left panel but with CO contours overlaid on the 
HST image obtained with NICMOS in the F160W band (color scale as shown in arbitrary units). Right panel Same as middle panel but zooming 
in on the continuum emission image of IRAS 17208-0014 at 221.1 GHz, in contours (1%, 2%, 5%, 10%, to 90% in steps of 10% of the peak 
value = 30.4 mJy beam“^ - 250cr), overlaid on the HST NICMOS image. The filled ellipses represent the CO beam size. 


4.6. IRAS 17208-0014: kinematic modeling 

4.6.1 . Kinemetry fit: the coplanar solution 

As for NGC 1614, we used the CO mean-velocity field of 
IRAS 17208-0014 to find the best-fit Fourier decomposition 
of the line-of-sight velocities using kinemetry (Krajnovic et 
al. 2006). Figure 14 shows the radial profiles of the main terms 
of the Fourier decomposition, and Fig. 15 shows the rotating 
disk solution of kinemetry and the residuals of the fit. As ex¬ 
pected, circular rotation accounted for by the ci term is dominant 
throughout the disk (Fig. 14a). The best-fit geometry for the ro¬ 
tating disk is defined by a position angle PA = 113 + 3° and 
an inclination i = 55 + 3°. We also determined L>sys(HEL) = 
12808 + 15 km s“^ as part of the fit. Relative to ^circ, non¬ 


circular motions (fnonc) are within a 2% to 50% range with a 
clear trend that shows higher non-circular motions at larger radii 
(Fig. 14b). Furthermore, based on the sign of the term, the 
pattern of the purely radial motions is undefined in the central 
r ~ 172 (1 kpc) region, where the 5 ' 1 /cl ratio shows an oscil¬ 
lating profile compatible with zero. Outside this central region, 
motions are indicative of a moderate outflow in the radial range 
r ~ 172 - 179 (1-1.5 kpc), while they suggest inflow farther out 
at r > 179 (1.5 kpc). Deprojected, these motions would trans¬ 
late into a maximum amplitude of < 25 km s“7 The ^-i term in 
IRAS 17208-0014 is nevertheless a small fraction (< 0.2 - 0.3) 
of Cnonc in tho outer molecular disk (r > 175 (0.8 kpc)). 
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IRAS17208-0014 





Fig. 14. Same as Fig. 6, but derived from the Fourier decomposition of the velocity held of IRAS 17208-0014. 

IRAS17208-0014 



Fig. 15. Same as Fig. 7, but showing here the comparison between the CO(2-l) observations and the best-fit model in IRAS 17208-0014. Contours 
in the left panel span the range -200 km s~^ to 200 km s"^ in steps of 50 km s"^ The CO(2-l) intensity contours in the right panel are as in the 
left-hand panel of Fig. 12. Color scales in both panels are as shown. 


4.6.2. A non-coplanar outflow solution 

Overall, similar to NGC 1614, the order of magnitude of the 
fitted noncircular motions can be explained as due to density 
waves. However, this scenario, which implicitly assumes that 
gas is in a coplanar geometry, is unable to explain the order of 
magnitude or even the asymmetric pattern of the high velocities 
ascribed to the line wing, as argued below. 

The left-hand panel of Fig. 16 shows the major axis p-v plot 
(along PA = 113°) with the projected best-fit rotation curve su¬ 
perposed. A sizeable fraction of the emission in the line core lies 
within the virial range determined by a combination of circular 
rotation (~ +230 km s“^), turbulence (FWHM/2~ ±85 km s“^), 
and an upper limit to the likely contribution from in-plane non¬ 
circular motions^. While the left-hand panel of Fig. 16 shows 
that the model, limited by beam smearing effects, may have 
underestimated the rotation curve gradient in the central r < 
075 (0.4 kpc) region, a hypothetically steeper rotation curve 
would nevertheless fail to explain the line wing feature: the stark 
asymmetry of the line wing, only detected at blue velocities, ar¬ 
gues against rotational support as an explanation of high veloc- 

This amounts to +115 km s"^ in IRAS 17208-0014. 


ities. Furthermore, as illustrated best in the right-hand panel of 
Fig. 16, with the achieved spatial resolution we locate the line 
wing feature in a quadrant of the p-v plot that is formally as¬ 
signed to gas in (coplanar) counter-rotation. Nevertheless, the 
coexistence of two opposite spin gas components at the same 
radii is highly unlikely because of the intrinsic dissipative nature 
of the gas. 


The alternative non-coplanar solution for the line wing in 
IRAS 17208-0014 can be used to invoke either inflow or outflow, 
depending on the adopted geometry. However, the detection by 
Arribas et al (2014) of a similar, yet less extreme, blueshifted 
Ha component in IRAS 17208-0014 favors the outflow scenario 
for the line wing, similar to the case of NGC 1614 discussed in 
Sect. 4.3. In addition, the atomic and molecular outflows seen, 
respectively, in the Nal D line (Rupke & Veilleux 2013) and in 
the OH line (Sturm et al. 2011) give further support to the out¬ 
flow scenario as an explanation for the CO extreme velocities in 
IRAS 17208-0014. 
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Fig. 16. Left panel. The (p-v) plot taken along the kinematic major axis {PA = 113°) of IRAS 17208-0014. As in Fig. 8, we delineate the 
(projected) rotation curve and the allowed virial range around it. Similarly, we delimit the region where emission is detected from the line wing. 
Middle panel: The p-v plot taken along the kinematic minor axis {PA - 23°). As in Fig. 8, the allowed virial range is delimited by the blue curves. 
Right panel: Same as left panel but zooming in on the central +1" region around the center. In all panels, contour levels are -2cr, 2cr, 2.5cr, 3cr, 
5cr, lOcr, 20cr, to 220cr in steps of 20cr with Icr = 1.2 mJy beam"^ The velocity scale is relative to v^y^. The spatial scales (X and Y) are relative 
to the dynamical center. We highlight the spatial resolution (beam size) achieved with UN weighting in each panel. 


5. Basic properties of the moiecuiar outfiows of 
NGC 1614 and IRAS 17208-0014 

5.1. Morphology 

Figure 17 shows the CO intensity maps of the molecular out¬ 
flows in NGC 1614 and IRAS 17208-0014, derived by integrat¬ 
ing the emission of the line wing components defined, respec¬ 
tively, in Sects. 4.1 and 4.4. We cannot rule out that emission at 
lower radial velocities apparently within the virial range is not 
coming also from outfiowing gas seen at a larger angle relative 
to the line of sight. However, we purposely adopt here the most 
conservative definition of the outfiows by restricting their veloc¬ 
ity ranges to the most extreme values associated with the line 
wings (see however discussion in Sect. 5.4). 

The molecular outfiow in NGC 1614 consists mainly of 
two barely resolved knots of emission located close to the nu¬ 
cleus along PA ~ 90° at {Aa, A^)=(-h0'.'8 ~ 250 pc, 0") and 

(-175 -470 pc, 0") (see the left-hand panel of Fig. 17). The 

east (west) knot corresponds to the blueshifted (redshifted) lobe 
of the outflow. The knots are bridged by lower level emission at 
the center. This morphology is reminiscent of a biconical out¬ 
flow uniformly filled by the outflowing gas. Overall, the bulk of 
the molecular outflow is cospatial with the 2"-diameter Patr star¬ 
forming ring imaged by Alonso-Herrero et al. (2001). However, 
the Ha emission coming from the blueshifted lobe of the out¬ 
flow, identified by Bellocchi et al. (2012), extends farther out up 
to r ~ 5" (1.6 kpc) on the eastern side of the disk. In this re¬ 
gion, there is tentative evidence of outflowing CO emission that 
stems from a number (~ 4) of clumps. The velocity range of the 
outflow is comparable in CO and Ha: the approaching outflow 
lobe is blueshifted by ~ 200 - 500 km s“^ relative to ^sys in both 
tracers. 

The molecular outflow in IRAS 17208-0014, identified as a 
blueshifted component, is associated with an off-centered knot 
of barely resolved emission located at (Aa, A^)=(0, -072 ~ 
-160 pc) (see the right-hand panel of Fig. 17). Within the er¬ 
rors this coincides with the position of the peak H 2 emission 
mapped by Medling et al. (2014), attributed to shocked gas at 


the base of an outflow coming from the nucleus of one of the 
two stellar disks (Medling et al. 2015). In contrast, the neutral 
atomic outflow traced by the blueshifted Nal D absorption comes 
from a wide angle (> 50°) and spatially extended component 
(up to r > 2" ~ 1.6 kpc in the plane of the sky) (Rupke & 
Veilleux 2013). Overall, the morphology of the atomic outflow 
is elongated along the CO kinematic major axis with hardly any 
extension along the minor axis, mimicking a wide angle (single) 
conical outflow. The CO outflow shows radial velocities simi¬ 
lar to those measured for the atomic gas in the inner r ~ 1" 
(0.8 kpc): velocities are blueshifted by ~ 600 km s“^ relative to 
fsys in both tracers. However, the (blueshifted) velocities of the 
H 2 outflow measured by Medling et al. (2015) are ~a factor of 
5 lower than the CO values. In Sect. 5.4 we propose a redefini¬ 
tion of the line wing velocities, which results in a change in the 
morphology, mass, and energetics of the CO outflow. 

5.2. Mass, energy, and momentum rates 

To quantify the mass load budget (Mout we use the 

expression 

^out — 3 X Vout ^ ^^outZ-^out ^ tan(a), (1) 

where Mout is the molecular gas mass of the outflow, Rout is its 
(projected) radial size, and Vout is the (projected) velocity of the 
outflowing gas. We also have to assume a certain geometry given 
by the angle a, which stands for the angle between the outflow 
and the line of sight. Equation (1) assumes a conical outflow uni¬ 
formly filled by the outflowing gas (Maiolino et al. 2012; Cicone 
et al. 2014). 

Alternatively, in the case of a single explosive (not continu¬ 
ous) event, an estimate of the outflow rate can be derived from 

^out — Eout ^ ^outZ-^out ^ tan(a), (2) 

where the ratio (VoutZ^out)~^ represents the dynamical time re¬ 
quired for the gas to reach its present location. The outflow rate 
from Eq. (2) is one third of the value derived from Eq. (1). 
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Fig. 17. Left panel: Overlay of the CO(l-O) intensity contours of the disk of NGC 1614, obtained by integrating the emission of the NA weighting 
data set inside the line core, on the intensity map of the molecular outflow (in color), obtained by integrating the emission inside the line wing 
region. The two velocity intervals are deflned in Sect. 4.1. Color scale spans the range [5cr,llcr], and contour levels are: 5%, 10%, 15%, 20%, 
to 90% in steps of 10% of the peak value = 38 Jy km s~^ beam"^ Right panel: Same as left panel but obtained from the UN weighting data set 
of IRAS 17208-0014. The color scale spans the range [5cr,15cr], and contour levels are 1%, 2%, 5%, 10%, to 90% in steps of 10% of the peak 
value = 127 Jy km s"^ beam"C The two velocity intervals are deflned in Sect. 4.4. 


Table 1. Molecular outflow properties. 


Sources 

^out 

Vout 

^out 

^out 

^out 

SFR 

^out 

Fut 

^SNe 

Fut 

r^AGN 

^out 

P out 
Lhol / C 

Pout 

Lagn/c 


(pc) 

(km s”^) 

(Me) 

(Me yr-*) 

- 

(Le) 

- 

- 

(g cm s"’) 

- 

- 

NGC 1614 

560 

360 

3.2 X 10’ 

40 

0.8 

6.5 X 10* 

0.06 

> 0.05 ^ 

1.4 X 10*5 

1.9 

> 65“ 

IRAS 17208 

160 

600 

4.6 X 10’ 

330 

1.4* 

1.6 X 10^0 

0.4* 

>6“^ 

1.6 X 10*5 

5 

> 5 X 10*" 


- 

- 

- 

- 

2.0 

- 

0.6 

0.02" 

- 

- 

17" 


Notes. All properties are derived for the line wing components as deflned in Sects. 4.1 and 4.4 assuming a - 33° and 35° in Eqs. (1), (3), and (4) 
for NGC 1614 and IRAS 17208-0014, respectively. Outflow mass rates, momentum rates, and kinetic luminosities have uncertainties of +0.34 dex, 
+0.36 dex, and +0.40 dex, respectively. The mass, momentum, and energy loading factors have associated uncertainties of +0.43 dex. See Sect. 5.3 
for details. 

(a) Derived assuming < 1.3 x 10^^ L© (Pereira-Santaella et al. in prep.). 

Derived assuming 5 FR = 240 M© yr“^ 

Derived assuming < 2.4 X 10^ L© (Gonzalez-Martm et al. 2009). 

{d) Derived assuming 5FR = 0.7 x 240 = 168 M© yr"^ 

Derived assuming ~ 7.2 X 10^^ L© (Aalto et al. 2015b.). 


Equation (1) was used by Cicone et al. (2014) to study the 
properties of the molecular outflows in their compiled sample of 
galaxies. In the case of NGC 1614, the observed morphology of 
the outflow favors the biconical fllled outflow scenario, and we 
also adopt here Eq. (1) in IRAS 17208-0014 for consistency and 
to facilitate the comparison with the work of Cicone et al. (2014). 

The kinetic luminosity (Fout) and momentum rate 
(Fout =dFout/d0 of the outflows are derived, respectively, 
from Eqs. (3) and (4), as follows 

^out — 1/2 X X (Vout/ cos((r)) (3) 


^out ~ ^out ^ ^out/ cos(cr). (4) 

Table 1 lists the values of Mout , Vout, and Fout used to cal¬ 
culate Mout, ^out, and Fout in NGC 1614 and IRAS 17208-0014. 


The masses Mmoi, which include the mass of helium, were cal¬ 
culated from the CO data cubes by integrating the emission of 
the lines in the line wings and assuming a conservatively low 
CO-to-H 2 conversion factor for the CO(l-O) line, which is typ¬ 
ical of mergers: ~ 1/5 of the MW value; i.e., we adopt Xqo = 
1/5 X 2 X 10^^cm“^(K km s“^)“^ (Downes & Solomon 1998; 
Tacconi et al. 2008; Daddi et al. 2010; Genzel et al. 2010; Bo- 
latto et al. 2013a). Eor IRAS 17208-0014, we adopt a 2-1/1-0 
brightness temperature ratio of 1 to derive Mout front the CO(2- 
1) integrated flux. 

The implied outflow rates given by Eq. (1) are Mout ~ 
62 X tan(Qf) M© yr“^ in NGC 1614, and ~ 510 x tan(Q') M© yr“^ 
in IRAS 17208-0014. Assuming that the outflow in NGC 1614 
is perpendicular to the main disk, a geometry suggested by the 
observed kinematics discussed in Sect. 4.1, cr = i = 33°, then 
^out ~ 40 M© yr“^ We have discarded the coplanar solution 
for the IRAS 17208-0014 outflow {a = 90° - i = 35°, with 
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i = 55°) in Sect. 4.6, but its geometry is less well constrained. 
The observed kinematics nevertheless suggest that the approach¬ 
ing cone is not perpendicular to the disk (a < i = 55°). We there¬ 
fore adopt 35° in the following as an indicative upper limit to a 
in IRAS 17208-0014, which implies that Mout ~ 330 M© yr“^. 

5.3. The powering sources of the molecular outflows: SF or 
AGN -driven? 

To investigate what the likely drivers of the observed molecular 
outflows are (star formation and/or AGN feedback), we derived 
the mass, energy, and momentum loading factors. These are de- 
flned, respectively, as: 1) the ratio of Mout to the integrated S FR 
(Mout/ S FR), 2) the ratio of Lout to the estimated luminosity of 
the AGN (Lout/LAGN) and to the total kinetic luminosity injected 
by supernova explosions (Lout/LsNe), and 3) the ratio of Lout to 
the total momentum rate transferred either by the AGN photons 
(^out/(^AGN/c)) or by all the photons (Pout/(Lho\/c)). The kinetic 
luminosity due to supemovae is estimated as LsNe(org s“^) ~ 7 x 
10^1 5 FR (Moyr-i)(Veilleux et al. 2005). We evaluate Lboi ~ Ljr 
in NGC 1614 and Lboi ~ 1.15 x Ljr in IRAS 17208-0014, in line 
with the prescriptions adopted by Veilleux et al. (2009) and Ci- 
cone et al. (2014) for luminous and ultraluminous starbursts. 

We have explored a range of values for Lagn- In NGC 1614 
Herrero-Illana et al. (2014) estimated an upper limit to the non- 
detected AGN in this source LAGN/^boi < ^0%, i.e., Lagn < 
4.5 X 10^^ L©. In the following we nevertheless adopt the more 
restrictive upper limit to the AGN luminosity in NGC 1614 de¬ 
rived by Pereira-Santaella et al. (in prep.), who combined the 
continuum observations at 435yum of ALMA (Xu et al. 2015) 
and the high-resolution MIR observations (at 24.5yum) done with 
the CanariCam instrument of the GTC, to do a SED fitting of 
the nuclear emission in a O'.'5-aperture: Lagn < 1.3 x lO^^L© ~ 
3%Lboi. In IRAS 17208-0014, the values discussed in the lit¬ 
erature for Lagn span a wide range from LAGN/^boi < 0.01%, 
i.e., Lagn < 2.4 x 10^ L©, as determined by Gonzalez-Martm et 
al. (2009), based on X-ray observations, up to LAGN/^boi ~ 0-3; 
i.e., Lagn ~ 7.2 x lO^^L©, as derived by Aalto et al. (2015b), 
based on the modeling of the vibrationally excited HCN emis¬ 
sion detected in the nucleus. Table 1 lists the values obtained for 
the loading factors as defined above. 

We have evaluated the uncertainties associated with the 
quantities listed in Table 1 (and also those listed in Table 2 and 
the related Fig. 20). Outflow mass rates, momentum rates, and 
kinetic luminosities have estimated uncertainties of +0.34 dex, 
±0.36 dex, and +0.40 dex, respectively. These include the 
+0.3 dex uncertainties on the mass estimates caused by calibra¬ 
tion and conversion factors, as well as the errors on the sizes, 
the inclinations, and the velocities of the outflow, contributing 
+0.1 dex each. The corresponding mass, momentum, and en¬ 
ergy loading factors have associated uncertainties of +0.43 dex, 
assuming that uncertainties on S FR and Lsnc. derived from Ljr, 
are of +0.17 dex. The errors on Lagn (and those derived for 
the magnitudes that depend on them) are reflected in the up¬ 
per limit and wide range of values available for NGC 1614 and 
IRAS 17208-0014, respectively. Overall, these error estimates 
are in line with the ones commonly derived in the literature for 
observations of molecular outflows similar to those presented in 
this work (see, e.g., Cicone et al. 2014). 

An inspection of the results listed in Table 1 leads to the fol¬ 
lowing conclusions regarding the nature of the powering sources 
of the outflows in NGC 1614 and IRAS 17208-0014. 

—NGC 1614: 


The bulk of the total SFR comes from the nuclear ring 
(~ 50 M©yr“^; Alonso-Herrero et al. 2001; U et al. 2012). In this 
region, where the molecular outflow is detected, S FR is compa¬ 
rable to Mout- Furthermore, the energy and momentum require¬ 
ments of the outflow, measured by Lout/LsNe and Pout/(Lho\/c), 
can be met by the star formation activity of the ring, assum¬ 
ing a conservatively low coupling efficiency (~ 6%) for the 
energy injected by supemovae explosions and a very low mo¬ 
mentum boost (~ 1.9) for all the photons (which come mainly 
from star formation). Similarly, based on the most restrictive 
limit on Lagn, the putative AGN in NGC 1614 could also drive 
the outflow, although in this scenario the momentum boost of 
AGN photons should be > 60. This factor is nevertheless at 
the higher end of the range of values predicted by AGN feed¬ 
back models under the assumption that molecular outflows are 
energy-conserving: Lout/(7^AGN/<^) ~ 10-50 (Faucher-Giguere 
& Quataert 2012). In summary, both star formation and AGN 
activity could cooperate to drive the molecular outflow in this 
source. The geometry of the line wing emission described in 
Sect. 4.1 indicates that the outflow axis is nearly perpendicular 
to the rotating disk, and it is spatially extended. While we cannot 
exclude that an AGN can also contribute to launch the outflow in 
a direction perpendicular to the disk, that the region in the disk 
where the outflow is detected extends on scales ~ 0.5 - 1 kpc 
suggests that star formation might be its main driving agent. 

—IRAS 17208-0014: 

As in NGC 1614, the integrated S FR ~ 240 M©yr“^ derived 
from a star-formation-dominated IR luminosity using a Chabrier 
IMF and the prescription of Kennicutt (1998), is comparable to 
Mout- Furthermore, the energy requirements of the outflow can 
only be met by assuming a significantly high coupling efficiency 
(~ 40%) for the energy injected by supemovae. The momentum 
boost factor is nevertheless moderate (~ 5) if all photons are con¬ 
sidered to contribute. In this scenario where most of the bolomet- 
ric luminosity comes from star formation (Lagn < 2.4 x 10^ L©; 
Gonzalez-Martm et al. 2009), the AGN power falls short of ac¬ 
counting for the luminosity and momentum rate of the molecu¬ 
lar outflow by large factors. However, if the limit on the AGN 
luminosity is as high as ~ 30% of Lboi (Aalto et al. 2015b), 
i.e., a factor of three higher than estimated from MIR diagnos¬ 
tics by Rupke & Veilleux (2013), this hidden AGN would be 
able to explain the molecular outflow with a low coupling ef¬ 
ficiency, Loui/Lagn ~ 0-02, and a moderate momentum boost 
factor. Pout/ (Lagn /r^) ~ 17. The geometry of the fine wing emis¬ 
sion discussed in Sect. 4.4, which indicates that the outflow axis 
is not perpendicular to the large-scale rotating structure, suggests 
that the outflow is launched by a non-coplanar disk that hosts ei¬ 
ther a compact starburst and/or an AGN. A merger episode could 
explain a random orientation of the inner circumnuclear gas disk 
relative to the larger scale disk. We describe in Sect. 5.4 an alter¬ 
native (less restrictive) definition of the velocities of the outflow 
in IRAS 17208-0014 and the implications thereof. 


5.4. An extreme AGN-driven outflow in IRAS 17208-0014? 

The molecular outflow axis in IRAS 17208-0014 is significantly 
tilted relative to the kinematic major axis of the rotating disk, 
as discussed in Sect. 4.6.1. To evaluate the intrinsic velocity 
boost in the outflow relative to the rotating frame of the disk in 
IRAS 17208-0014, we obtained a new version of the data cube 
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Fig. 18. CO(2-l) velocity-channel maps of IRAS 17208-0014 obtained after subtracting the rotation curve model, as discussed in Sect. 5.4. 
Contours, symbols, and scales are the same as in Fig. 9. The dividing line between the line core and the line wing is here shifted to \v - i^sysl ~ 
300 km s“^ 


IRAS17208-0014 



in 




Fig. 19. Same as Fig. 11 but for the centroids of CO(2-l) emission in IRAS 17208-0014 derived from the data cube after subtraction of the rotation 
curve. Red is used for velocity channels of the line wings as redeflned in Sect. 5.4. Black is used for the line core. 


by subtracting the rotation curve model derived in Sect. 4.6.1^ at 
each pixel. 

Figure 18 shows the same velocity-channel maps as are dis¬ 
played in Fig. 9 but after subtracting the rotation curve. This 


^ We note that applying the same correction to the data cube of 
NGC 1614 (results are not shown here) had no signiflcant effect on the 
velocity pattern of the outflow because the latter is oriented close to the 
minor axis of the rotating disk, as discussed in Sect. 4.3.1. 


figure highlights that the decoupling of the outflow component 
relative to the rotating disk, which is revealed by the apparent 
reversal of the respective kinematic major axes, is identified by 
visual inspection at | - t>sys I ^ 300 km s“^ This result signifi¬ 

cantly expands the velocity range of the line wings, here defined 
relative to the rotating frame of the disk, compared to the de¬ 
scription of Sect. 4.4. In particular, the redshifted counterpart of 
the outflow, which remained unnoticed in Fig. 9, is now identi¬ 
fied in the range v - v^ys = [+300, +540] km s“^ (see Fig. 18). 
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Table 2. Re-evaluation of the moleeular outflow properties in IRAS 17208-0014 


Source 

Fout 

Vout 

3Fout 

3Fout 


T 

-^out 

-^out 

P 

Pout 

Pout 

SFR 

^OUt 

^SNe 

r^AGN 

^out 

^bol / C 

Lagn/c 


(pc) 

(km s“^) 

(Mo) 

(Mo yr-‘) 

- 

(A) 

- 

- 

(g cm s"^) 

- 

- 

IRAS 17208 

160 

500 

2x 10^ 

1200 

5^ 

4 X 10*0 

1 ^ 

> 15*' 

4.8 X 10^0 

15 

> 1.5 X 10^* 


- 

- 

- 

- 

7" 

- 

1.4^ 

0.05 

- 

- 

51 ‘‘ 


Notes. New properties derived for the line wings as deflned in Sect. 5.4 assuming a = 35°. Outflow mass rates, momentum rates, and kinetic 
luminosities have uncertainties of +0.34 dex, +0.36 dex, and +0.40 dex, respectively. The mass, momentum, and energy loading factors have 
associated uncertainties of +0.43 dex. See Seet. 5.3 for details. 

(a) assuming S FR = 240 M© yr“^ 

Derived assuming < 2.4 x 10^ L© (Gonzalez-Martm et al. 2009). 

Derived assuming 5 FR = 0.7 X 240 = 168 M© yr“^ 

(d) Derived assuming ~ 7.2 X 10^^ L© (Aalto et al. 2015b). 


To quantify the effect of the decoupling, we show the cen¬ 
troids of CO emission in IRAS 17208-0014 in Fig. 19, as well as 
the PA and r profiles obtained for the velocity channels shown in 
Fig. 18. As illustrated by the middle panel of Fig. 19, the value 
of PA for I i; - t>sys I ^ 200 km s“^ is FAcore = 107 + 8°, i.e., close 
to the PA of the rotating disk derived in Sect. 4.6.1. This indi¬ 
cates that the rotation pattern is still present in this velocity range 
of the kinematically rescaled data, reflecting that the model un¬ 
derestimates the intrinsic rotation curve, as already mentioned 
in Sect. 4.6.2. In stark contrast, the PA values for the centroids 
at I f - L>sys I ^ 300 km s“^ oscillate around a weighted mean of 
FAwing = 308° +5°, revealing a true reversal of the kinematic ma¬ 
jor axis of the outflow relative to the disk. We conclude that the 
evidence of a kinematic decoupling, reflected in the measured 
difference FA wing-FAcore ~ 201+9° is a significant > 20(T result. 
Furthermore, the line wing channels show slightly lower r values 
(~ 072 (160 pc)) compared to the line core defined in Sect. 4.4 
(up to ~ 073 (240 pc); see the right-hand panel of Fig. 11). 

Table 2 lists the new values obtained for the loading factors 
of the outflow in IRAS 17208-0014 derived from the redefined 
velocity range. The enhanced mass loading factor (Mout/5' FR ~ 
5-7) cannot be easily explained by star formation activity 
alone. From the observational point of view, mass loading factors 
exceeding unity are rarely seen in actively star-forming merg¬ 
ers, in particular in the range SFR > 100 M© yr“^ applicable 
to IRAS 17208-0014 (Rupke et al. 2005b; Sturm et al. 2011). 
Moreover, Table 2 shows that the required efficiency for convert¬ 
ing the energy injected by supemovae into mechanical energy 
of the outflow would be unrealistically high (> 100%). Charac¬ 
teristic values defined observationally for the energy conversion 
efficiency oscillate between 10% and, at most, 50% (Veilleux et 
al. 2005; Veilleux 2008). By contrast, a putative hidden AGN of 
Fagn ~ 7.2 X 10^^ L© would be a viable driving agent of the 
molecular outflow in this source. This scenario would require a 
moderate coupling efficiency, Lqui/Fagn ~ 0-05, and an admit¬ 
tedly high momentum boost factor, Foui/CFagn/^^) ~ 51. These 
values are nevertheless in line with those commonly derived in 
other AGN-driven outflows (Cicone et al. 2014). 

Independent evidence of a possible hidden AGN in 
IRAS 17208-0014 can be found in the new ALMA obser¬ 
vations of the {v2 = 1) 4-3 HCN vibrational line emission 
(Aalto et al. 2015b). IRAS 17208-0014 is part of a family 
of (U)LIRGs where the emission in vibrationally excited lines 
seems to emerge from buried, compact (r < 17-70 pc) nu¬ 
clei that have very high implied mid-infrared surface brightness 
> 5 X 10^^L© kpc“^. These nuclei are likely powered by accret¬ 


ing supermassive black holes (SMBHs) and/or hot (> 200 K) 
extreme starbursts. The reported extreme surface brightness of 
the hidden core in IRAS 17208-0014 would be able to excite the 
HCN vibrational line by intense 14 jim emission. The new obser¬ 
vations of IRAS 17208-0014 presented by Aalto et al. (2015b) 
also show evidence of an extreme outflow in the dense molecular 
gas traced by the ground vibrational state {v = 0) HCN(4-3) line 
emission imaged by ALMA. 

6. The molecular outflows of NGC 1614 and 
IRAS 17208-0014 in context 

We compiled a sample of 23 nearby galaxies where molecular 
outflows have been observed and detected to date, taking the 
sample of 18 objects studied by Cicone et al. (2014) as start¬ 
ing point, complemented here by new data recently obtained in 
five targets. These additions include two Seyferts (NGC 1433: 
Combes et al. 2013 and the CND of NGC 1068 = NGC 1068- 
CND: Garcfa-Burillo et al. 2014) and three mergers (the southern 
nucleus of NGC 3256 = NGC 3256-S: Sakamoto et al. 2014 and 
Emonts et al. 2014, as well as NGC 1614 and IRAS 17208-0014: 
this work). We have also included the new outflow parameters 
derived by Alatalo et al. (2015) in NGC 1266. Figure 20 illus¬ 
trates the different scaling relations (on logarithmic scale) ob¬ 
tained in an attempt to identify the driving agents of the molec¬ 
ular outflows inside the parameter space probed by the sample. 

The new data of NGC 1433, NGC 1068-CND, NGC 1614, 
and NGC 3256-S follow the correlation between the outflow 
rate and the AGN luminosity found by Cicone et al. (2014) in 
their sample, which is an indication that the molecular outflow 
in these sources can be AGN-driven (left-hand panel of Fig. 20). 
The location of IRAS 17208-0014 in the SFR - Fagn diagram 
depends critically on the wide range of values for the AGN 
power discussed in Sect. 5.3. This source follows the correla¬ 
tion if the highest value of Fagn is adopted but it becomes an 
outlier, like the pure starbursts of Cicone et al. (2014), if the low 
end of Fagn is assumed. 

The middle panel of Fig. 20 shows the outflow rate as a func¬ 
tion of the SFR. In this diagram, NGC 1433, NGC 1068-CND, 
and NGC 3256-S share the location of other AGN-driven out¬ 
flows as they all depart from the 1:1 line followed by pure star- 
bursts. In contrast, NGC 1614 follows the correlation of pure 
starbursts, an indication that star formation can also cooperate 
to drive the outflow in this source. The location of IRAS 17208- 
0014 in this diagram depends on the value adopted for Mout- the 
more extreme loading factors discussed in Sect. 5.4 put it be- 
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Fig. 20. Scaling relations derived for the molecular outflows observed in the sample of galaxies compiled for this work: outflow mass-loss rate 
as a function of the AGN bolometric luminosity (left panel), outflow mass-loss rate as a function of the SFR {middle panel), and momentum 
boost factor as a function of mass-loading factor {right panel). Open symbols identify galaxies from the sample of Cicone et al. (2014). Open 
(red) triangles show the data from pure starburst galaxies; open (black) squares identify the data of low and high-luminosity AGNs, where the 
presence of an AGN has been securely detected. Filled symbols identify the data from the new additions discussed in Sect 6: NGC 1433 (Combes 
et al. 2013), the CND of NGC 1068 = NGC 1068-CND (Garcia-Burillo et al. 2014), the southern nucleus of NGC 3256 = NGC 3256-S (Sakamoto 
et al. 2014; Emonts et al. 2014), as well as NGC 1614 and IRAS 17208-0014 (this work). The hashed polygons show the allowed range of values 
for IRAS 17208-0014 in all the panels. The line in the left panel shows the linear fit to the AGN sample. Lines in the middle and right panels 
illustrate different values of the mass and momentum loading factors. One-sided errorbars identify upper and lower limits on the values plotted in 
all the panels. In the bottom right of the left and middle panels and in the bottom left of the right panel we show Icr errorbars that are representative 
of the whole sample used in this work (see discussion in Sect. 5.3 and Cicone et al. 2014). 


yond the 2-3:1 line predicted by numerical simulations that in¬ 
clude radiative and mechanical feedback in star-forming galaxies 
(Hopkins et al. 2012; Lagos et al. 2013). 

The right-hand panel of Fig. 20 shows the momentum boost 
factor as a function of the mass-loading factor for the compiled 
sample. Galaxy outflows in the upper left-hand (lower right- 
hand) quadrant with momentum-boost factors > 50 (< 50) and 
mass-loading factors < 2.5 (> 2.5) qualify as star formation 
(AGN) driven systems. Galaxies lying in the lower left-hand 
quadrant are mixed systems to the extent that both star for¬ 
mation and AGN can cooperate to drive the outflow in these 
sources. NGC 1614 is close to this undefined region. Galaxies 
like NGC 3256-S, which lies in the upper right-hand quadrant, 
likely require other mechanisms like radio jets to drive the out¬ 
flow. The location of IRAS 17208-0014 in this diagnostic dia¬ 
gram reflects a range of possible classifications. However, if we 
assume the properties derived in Sect. 5.4, a hidden AGN is the 
most likely powering source of the outflow in this source. 

7. Summary and conclusions 

We have used the PdBI to image the CO emission in two- 
starbursts dominated (U)LIRGs: NGC 1614 and IRAS 17208- 
0014 with high spatial resolution (075-172). The two targets, 
classified as mergers in an advanced stage of interaction, are 
vigorous star-forming systems {SFR > 50 - 240 M© yr“^), and 
they present evidence of an outflow in several ISM phases, in¬ 
cluding the ionized and the neutral gas. The high resolution of 
the PdBI observations used in this work are instrumental in con¬ 
straining the geometry, as well as the mass, momentum, and en¬ 
ergy budgets of the molecular outflows discovered in NGC 1614 
and IRAS 17208-0014. 

We summarize the main results of our study as follows: 

- The CO(l-O) emission in NGC 1614 has been mapped with 

a spatial resolution of 175 x 078. The bulk of the CO emis¬ 


sion comes from a spatially resolved (3.4 kpc x 2.5 kpc), 
lopsided molecular disk. The disk shows an arc-like feature 
that runs from the northeastern to the southeastern side. A 
sizable fraction of the CO emission originates in the 600 pc- 
diameter star-forming ring of supergiant HII regions, which 
is associated with a young starburst {SFR ~ 50 M© yr“^). 
Fainter CO emission is also detected and extends farther out 
to the east up to r ~ 1.9 kpc. 

- The CO emission in NGC 1614 extends over a symmetric 
~ 900 km s“^-wide velocity span. The bulk of the emission 
comes from a more restricted velocity range: | v - v^ys I < 
210 km s“^ (the line core). We also detect emission that ex¬ 
tends up to I f - fsys I ~ 450 km s“^ (the line wings). Emis¬ 
sion from the line core shows the signature of rotating disk 
with a kinematic major axis oriented roughly north-south, as 
confirmed by a Fourier decomposition of the mean velocity 
field. In contrast, the kinematic major axis for the line wings 
is tilted progressively toward the east-west axis, an indica¬ 
tion of strong kinematic decoupling of the high-velocity gas. 
The kinematic pattern of the line wings can be explained by 
a non-coplanar molecular outflow. 

- The molecular outflow in NGC 1614 has a mass of 3.2 x 
lO^M© and consists mainly of two knots of emission that 
are co-spatial with the nuclear star-forming ring. The eastern 
(western) knot is blue(red)-shifted by ~360 km s“^ on aver¬ 
age. The blueshifted velocity ranges of the outflow are com¬ 
parable in CO and Ha. The mass, energy, and momentum 
budget requirements of the molecular outflow in NGC 1614 
can be met by the star formation activity taking place in the 
ring. The putative AGN in NGC 1614 could also drive the 
outflow. The geometry of the line wing emission indicates 
that the outflow axis is nearly perpendicular to the rotating 
disk, and it is spatially extended on scales ~ 0.5 - 1 kpc. 
This suggests that star formation is the most likely driving 
agent of the outflow, although a contribution from a putative 
AGN cannot be excluded. 
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- The CO(2-l) emission in IRAS 17208-0014 has been 
mapped with a spatial resolution 076 x 075. The CO emis¬ 
sion comes from a spatially resolved (2.7 kpc x 1.8 kpc) 
molecular disk. The morphology of the CO disk is simi¬ 
lar to the northeast-southwest elongated structure seen in 
the dust extinction, as derived from the 2.2/jm/l.lpm HST 
color images. Three fainter protrusions and a detached clump 
extend the disk emission farther out up to r ~ 1.8 kpc. 
The molecular disk feeds a vigorous star formation episode 
(SFR ~ 240 Mq yr-i). 

- The CO emission in IRAS 17208-0014 covers a wide and 
asymmetric velocity range: v - t>sys = [+450,-750] km s“^. 
However, up to ^ 99% of the total emission comes from a 
narrower velocity range: v - fsys ~ [+450, -450] km s“^ (the 
line core) that shows the characteristic feature of spatially re¬ 
solved rotating disk with PA = 113 + 3°, as determined by a 
Fourier decomposition of the mean velocity field. Rotation is 
perturbed by strong non-circular motions. Gas emission from 
the blueshifted line wing is kinematically decoupled from the 
disk: its apparent kinematic major axis is virtually reversed. 
This pattern can be explained by a non-coplanar molecular 
outfiow. This scenario is reinforced after subtracting the ro¬ 
tation curve, which reveals that the kinematic decoupling of 
the outfiow starts at | i; - t>sys I ^ 300 km s“^ 

- The molecular outfiow in IRAS 17208-0014 has a mass of 
2 X 10^ Mq and consists of two knots of emission located 
at r ~ 160 pc. The southeast (northwest) knot is blue(red)- 
shifted by -500 km s“^ on average. The blueshifted veloc¬ 
ity ranges of the outfiow are similar in atomic, ionized, and 
molecular gas. In stark contrast to NGC 1614, the mass, en¬ 
ergy, and momentum budget requirements of the molecular 
outfiow in IRAS 17208-0014 cannot be met by its star for¬ 
mation activity. Instead, a putative hidden AGN of Lagn ~ 
7.2 X 10^^ L© would be a viable driving agent of the molec¬ 
ular outfiow in this source. The geometry of the molecular 
outfiow, with an axis that is not perpendicular to the large- 
scale rotating structure, suggests that the outfiow is launched 
by a non-coplanar disk. 

The molecular outfiow in NGC 1614 is likely triggered by 
its vigorous star formation activity, which is mainly concen¬ 
trated in the 600 pc-diameter nuclear ring of the galaxy. This 
explains the extended nature of the outfiow, but also its orien¬ 
tation, which is roughly orthogonal relative to the large-scale 
rotating disk. In the case of IRAS 17208-0014, the molecular 
outfiow geometry and energetics can be explained best as driven 
by an AGN. It is tempting to identify which of the two over¬ 
lapping nuclear stellar disks discovered in IRAS 17208-0014, 
which lie 0'!2 (200 pc) apart and which were recently discov¬ 
ered by Medling et al. (2014), is responsible for launching the 
outfiow in this source. Although our spatial resolution (075) is 
insufficient to precisely locate the origin of the outfiow spatially, 
the observed kinematics can be used to make a tentative identifi¬ 
cation. Medling et al. (2014) characterize the two rotating stellar 
disks kinematically: the major axis is oriented at PA - 40° for 
the W nucleus and at PA - 128° for the E nucleus. While the 
PA of the large-scale rotating molecular disk (PA - 113°) is 
closer to the orientation of the E disk, the derived molecular out¬ 
fiow axis (PA - 308°) is, within a 2° tolerance, perpendicular to 
the W disk. This suggests that the launching mechanism of the 
molecular outfiow in IRAS 17208-0014 lies in a buried AGN 
associated with the W nucleus. To be obtained with the ALMA 
interferometer, higher resolution observations are required to un¬ 
ambiguously identify the launching mechanism of the molecular 
outfiow in IRAS 17208-0014. 
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